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Abstract A number of misunderstandings about modeling the apparent accelerated expansion of the Universe, 

I and about the 'weak singularity' are clarified: 1. Of the five definitions of the deceleration parameter given by 

■ Hirata and Seljak (HS), only qi is a correct invariant measure of acceleration/deceleration of expansion. The 
' (?3 and q4 are unrelated to acceleration in an inhomogeneous model. 2. The averaging over directions involved 

, in the definition of 54 does not correspond to what is done in observational astronomy. 3. HS's equation (38) 

^ ' connecting 54 to the flow invariants gives self-contradictory results when applied at the centre of symmetry of the 
Lemaitre-Tolman (L-T) model. The intermediate equation (31) that determines qy is correct, but approximate, 

^ ■ so it cannot be used for determining the sign of the deceleration parameter. Even so, at the centre of symmetry 

{^j^l of the L-T model, it puts no limitation on the sign of g3'(0). 4. The 'weak singularity' of Vanderveld et al. is 

' ' . a conical proflle of mass density at the centre - a perfectly acceptable conflguration. 5. The so-called 'critical 

' point' in the equations of the 'inverse problem' for a central observer in an L-T model is a manifestation of the 

^ . apparent horizon - a common property of the past light cones in zero-lambda L-T models, perfectly manageable 

' if the equations are correctly integrated. 

o ■ 
^ ■ 

• ' 1 Motivation 

m ; 

' Vanderveld et al. [1] (abbreviated as VFW) claimed that there was a contradiction between two results concerning 

, the putative accelerated expansion of the Universe. Their reasoning was, in brief, this (italics mark quotations): 

' On the one hand, Hirata and Seljak [2] (abbreviated as HS) claimed to have proved that in a perfect fluid 

. 5^ I cosmological model that is geodesic, rotation-free and obeys the strong energy condition p + 3p > 0, a certain 

■ generalisation of the deceleration parameter q must be non-negative. But on the other hand, Iguchi et al. [3] 
$H ' (hereafter INN) did obtain simulated acceleration in Lemaitre - Tolman (L-T) models [4,5] with A — that 

obey HS's conditions. This contradtctton is resolved by showing that L-T models that simulate accelerated 
expansion also contain a weak singularity, and in this case the derivation of HS breaks down. In addition to this, 
there are other singularities that tend to arise in L-T models, and VFW have failed to Gnd any singularity-free 
models that agree with observations. 

This leaves the impression that physically acceptable inhomogeneous models are unable to account for obser- 
vations. It is shown here that this reasoning is not correct. In brief, our theses are the following: 
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1. Of the five definitions of the deceleration parameter given by HS, only qi is a correct invariant measure 
of deceleration of fluid expansion. Their derives from a Taylor expansion of the luminosity distance-redshift 
relation, Dl{z), in powers of z, done separately for each direction. Their qi is based on the angular average of a 
Taylor expansion of z{D[^), and is conceptually different from what is done in observational practice. Although 
all these definitions reduce to the familiar q in FLRW models, they have different values and distinct meanings in 
an inhomogeneous model. In particular q^i and 94 are not measures of deceleration in an inhomogeneous model. 

2. HS's 54 does not correctly represent the observers' deceleration parameter q. When observers assume that 
our real Universe is in the FLRW class, q is the same in all directions. But when they consider inhomogeneous 
and anisotropic models, they have to measure the distance-redshift relation for each direction separately rather 
than average the measured results over directions, as this would mean destroying useful information. Thus, HS's 
qi that represents the result of such averaging does not correspond to the observations actually done in relation 
to this quantity. 

3. HS's reasoning is correct (in the approximate sense) up to their eq. (31), which allows one to calculate 93/ 
at the centre of symmetry in a Lemaitre - Tolman (L-T) model. But, at a spherical centre, the final result is the 
opposite to what they intended - no limitation on the sign of 54(0) follows (see our Sees. 4 and 6). 

Thus, there is no contradiction involved in a decelerating inhomogeneous model imitating observational rela- 
tions of an accelerating FLRW model. Even so, it is incorrect to use an approximate equation to decide whether 
some quantity is positive or negative. 

4. What VFW call a weak singularity is not a singulalrity. ^ 

5. It is also not true that other singularities invalidate the L-T models. The differential equations for the 
inverse problem - given observational data functions, calculate the LT model that would give them - are indeed 
singular at the apparent horizon. VFW refer to this location as a critical point, and to this phenomenon as a 
pathology. In truth, this is simply the reconvergence of the observer's past null cone towards the Big Bang, that 
is a consequence of the decelerating cosmic expansion when lambda is zero, long known in the FLRW CEise, e.g 
[7-9]. Though some authors (INN and VFW among them) were unable to propagate their solutions through 
the 'critical point', this difficulty was overcome by Lu and Hellaby [10] and in fact the AH relation was used to 
provide extra information by McClure and Hellaby [11,12]. 

Since these misunderstandings have propagated into the literature, it is essential to resolve these issues. 

But the most important point to be stressed is this: what one wants to reproduce, in connection with the 
distant type la supernovae, is not the accelerated expansion of the Universe, but the observed luminosity distance- 
redshift relation. The apparent dimming of supernovae (first reported by Riess et al. [13] and Perlmutter et al. 
[14]) stems from a comparison between observations and the Einstein-de Sitter model, the 'old' standard model 
of the Universe. The first proposed explanation, originating from the assumption that the Universe should be 
Friedmannian, is an accelerated expansion. However, this acceleration is not the phenomenon to be explained, 
but a component of a particular explanation - one that assumes homogeneity on the scale in question. 



2 The Lemaitre - Tolman (L - T) model 

Since in the following we often refer to the L-T model, we summarise here the basic facts about it in the case 
A = 0. For more extended expositions see Refs. [15,16]. The metric of the L-T model is: 

D 2 

ds^ = dt^ - ^^^^dr^ - R\t,r){M^+ sin' ^ dip''), (2.1) 

where E{r) is an arbitrary function, and R{t, r) is determined by the integral of the Einstein equations: 

r/ = 2E{r) + 2M{r)/R, (2.2) 

M(r) being another arbitrary function. Equation (2.2) has the same algebraic form as one of the well-known 
Friedmann equations, except that here it contains arbitrary functions of r in place of arbitrary constants. The 
Friedmann limit follows when E — — fcr^/2, M = Mqt^ and R = rS{t) where k, Mq and S{t) are the corresponding 
Friedmann constants and scale factor. The solution of (2.2) may be written as 

t-tB{r)= I (2.3) 

where tgir) is one more arbitrary function called the bang-time function; in the Friedmann limit it is constant. 
The + sign apiilics for an expanding region, — applies for a collapsing region. The mass density is 

^ Moreover, their definition of 'weak singularity' does not agree with the meaning of the term 'singularity' in relativity, 
and is not supported by anything in the Tipler paper [6] they cite. 
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2M r 

8.Gp=^. (2.4) 

The pressure is zero, and so the matter (dust) particles move on geodesies. 

The equations determining the L-T model are covariant with the transformations r — ^ r' = /(r). Therefore, 
we may use such a transformation to give one of the functions (M, ig) a handpicked form, provided the chosen 
function is monotonic in the range under investigation. 

An incoming radial null geodesic is given by the differential equation 

? = -^^^' (2-5) 
dr ^1 + 2E{r) 

and its solution is denoted t = t{r). The general solution determines the past null cone (PNC) of the observer 
situated at the centre of symmetry. We use the hat " and the subscript A to indicate evaluation on the PNC. 
Then the redshift along the PNC, z{r), is given by [17,16]: 



1 Az 
1 + « dr 

and the area and luminosity distances are 



Ra 



(2.6) 



Da = R, Dl = {1 + zfR. (2.7) 

Equations (2.5) and (2.6) can be solved numerically once M{r), E{r) and tgir) are specified. 

An L-T model may possibly have a curvature singularity at the centre of symmetry. The singularity will be 
absent when there is no point mass there, and the mass density is finite. This will happen when the conditions 
given below are obeyed (see Ref. [16] for the derivation). 

Let r = rc be the radial coordinate of the centre of symmetry, where R{t,rc) = for all t > tB{rc), and let 

def 

p{t,rc) = ct{t) < oo. We choose the r-coordinate so that rc = and 

M = Mor^, (2.8) 
where Mq = constant. Then, in a neighbourhood of r — 0, the functions behave as follows: 

i^ = /3(^)r + 7^l(r•), E = 'yr^+TZ2{r), ^B = T + 7^o(^), (2.9) 
7 and t being constants, where the symbols 7la{r) will denote quantities with the property 

Inn = 0. (2.10) 



3 The definitions of deceleration by HS and their physical meaning 

Hirata and Seljak [2] give in their paper five alternative definitions of the deceleration parameter, three of which 
we briefiy recall here for reference. In all cases it is assumed that the spacetime is filled with a perfect fluid 
obeying Einstein's equations with zero cosmological constant. HS assume that the fluid moves geodesically, but 
for the beginning we shall drop this assumption. 

Their first definition of the deceleration parameter is: 

qi - 1 - u^'Hi.^/Hi^ Hi m'';^/3, (3.1) 

where Hi is the (local) Hubble parameter and is the velocity field of the fluid. Then, using these definitions, 
the Raychaudhuri equation may be rewritten as follows: 

Hi\i = i [-uT;^ + a^-'a^, - ^^'•'io^,,] + ^(p + 3p), (3.2) 

where ii^ is the acceleration of the fluid, a^iiy and ujfiiy are the shear and rotation tensors, respectively, p is the 
mass density and p is the pressure. It follows that with w'* = = co^i, and p + 3p > 0, gi must be non-negative. 
This result applies in particular to the L-T models (where ii^ = = u)^v and p = 0) and to the FLRW models 
(where ti^ = = w^j/ = cr^i/). 
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The definition (3.1) is invariant, applies in all geometries, is a generalization of the deceleration parameter 
long used in the FLRW geometries, and measures the relative change of the expansion scalar along the flow lines 
of the perfect fluid. With 51 > we have w'*J7i;^ < —H\^ < 0, so the expansion slows down toward the future 
of the comoving observers. This is the case in all cosmological models that obey the assumptions listed above. 

The quantities 53' and H^^i refer to a Taylor expansion of the luminosity distance-redshift relation z{D]^): 

z = H^.Dl - \h^'^{1 - q3,)DL' + OiDL^), (3.3) 

where Hy and qy are functions of the angle of observation. The deflnitions of H4 and 54 follow by averaging 
(3.3) over the full solid angle 4-k at each fixed D^: 

<z>4^= H^Dl - \ha^{'^ - 1i)DL^ + 0{Dl^). (3.4) 
These equations have the same algebraic form as the corresponding equation in an FLRW model. 



4 Problems with HS's equation (38) 

By considering a bundle of null geodesies converging to a point in space, Hirata and Seljak [2] derived an equation 
relating the quantities in (3.4) to the invariants of flow of a perfect fluid flUing the spacetime. Flanagan's (hereafter 
F, [18]) equation is more general, and it reads, in HS's notation, 

H^^i = ^(p + 3p) + i + l^^pa'^^ - oj^iscj''^ - 2M«;a] . (4.1) 

(HS's result is the subcase = = Wa/j.) 

Equation (4.1) has the appearance of being exact and covariant. Unfortunately, it is neither. As we show 
below, its derivations, both by HS and by Flanagan, contain approximations (some not explicitly spelled out) 
and involve averaging over directions, which produces coordinate-dependent results at locations corresponding 
to coordinate singularities. Approximations may be good for some purposes, but not when one intends to show 
that some quantity is positive. 

Below we note instances where HS and F introduced approximations in their reasoning. 



4.1 HS's equation (26). 

We have no objections to the reasoning of HS up to their eq. (25). Having established the relation between the 
luminosity distance Dj^ and the redshift z (their (24)), and between the radiation amplitude and the expansion 
scalar On of the radiation field (their (25)), HS then copy a relation between dn {0 in HS's notation) and the 
afflne parameter v from Ref. [19] (HS's eq. (26)). This requires a comment. 

To obtain this relation one considers another family of rays, originating at the centre of a source (call it a 
star for brevity). The star is assumed spherical and having the radius R. One of the rays hits the observer. The 
values of the affine parameter along this ray are w = at the observer, Vs at the surface of the star and {vs + Avs) 
at the centre of the star. The relation is then 

= , ^, . , + 0[v - {vs + Ave)]. (4.2) 

V — (Vs + Avs) 

In Ref. [19] this is obtained as an approximate solution of the Raychaudhuri equation for null geodesies, which is 

fe^6l„,^ + 2 (ffn'^ - OJn^^ + = -RfyykPk'^, (4.3) 

where 9n = k^;ti, (Jn and Un are, respectively, the expansion, shear and rotation of the null congruence k'^ = 
dx^/dv and Rp-f is the Ricci tensor in spacetime. 

Equation (4.2) without the second term on the right is the exact solution of (4.3) in the case an = cun = = 
Rp-fk^k'^ . Thus the approximation involved in (4.2) is that the contributions to k'^9^~^ from shear, rotation and 
lensing by matter are negligible compared to dn^- These are additional assumptions about light propagation in 
spacetime [an and LOn negligible) and about the spacetime itself (Rp^k'^k'' negligible). 

However, HS's (38) is used by VFW as if it were exact. When an ~ — Rp^k^'k''' hold exactly, the Goldberg 
- Sachs theorem [20] applies, which says that the spacetime must be algebraically special, with k'^ being the 
double principal null congruence of the Weyl tensor. The assumption Rp-^k'^k'^ = 0, when imposed on the L-T 
spacetime, reduces it to the Scliwarzschild solution. Therefore, HS's equation (38), when treated as exact, involves 
strong assumptions that don't apply to the L-T model, or even the FLRW model. 
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4.2 The two errors that cancelled each other 



The two errors shown below have no meaning for the final result because (in the approximation in which HS 
worked) they cancelled out in the end, but they might be confusing. HS's equation (27) should actually read 

Av. = (4.4) 

i.e. it is the inverse of what they wrote. Similarly, their equation (28) should actually read 

(whore A is the radiation amplitude), i.e. again the inverse of their result. Since HS assumed \Avs\ <C \vs\, the 
above is approximately equal to 

^(■Us) Vs 1+z 

which results in HS's eqs. (29) - (31) being right again. Since we will later refer to their (31), we copy it here: 
z = -Kijn'n^DL - 2 (^Kijn'n^y {Dl^ + ^ (kij + ify|fen'= + AKi'^Kki) n'n^ {DLf + 0{Dl^). (4.7) 

The meaning of the symbols in (4.7) is as follows. Since HS work under the assumption that w = for the 
cosmic fluid, comoving and synchronous coordinates exist in which the metric gr^^ has the properties qqo — 1, 
goi = 0, where = t is the time coordinate and x^, i = 1,2,3, are the space coordinates. Then hij = —Qij is 
the metric of a space of constant t, and Kij is the second fundamental form of this space; in these coordinates 

def 

Kij = — (1/2) (d/dt) hij = — (l/2)ftjj. Then the geodesic null vector k" can be normalised at the observer so 
that = 1 and its space components k^ ri' form a 3-dimensional unit vector, hijU^nP = 1. 



4.3 The problem with averaging 

There is a problem with averaging products like n^nP and ri^n-'n^ at coordinate singularities in space, and the 
centre of symmetry in L-T is a singularity of the spherical coordinates. For example, a general unit vector in 
Euclidean space attached at a point O, when referred to Cartesian coordinates and parametrised by spherical 
angles {a, 13), has the components (n^,n^,n^) = (sin a cos /3, sin a sin ^, cos a). In this case, the averages'^ are 
<n'n''>4^ = (1/3)5'^ <n'n-'n''>4^ = 0, <n'n-'n^n'>4^ = (l/5)(5'^*''5'''^. The first and last of these are special 
cases (in Euclidean space) of the covariant equations 

<n'n^>4^ = {l/3)h'^, (4.8) 
<n''n^n''>4T, = 0, (4.9) 
<n'n^n'=n'>4^ = (l/5)/i(*^/i'=') = (1/15) (h'^ h''^ + h'^h^'' + /i^'/i^'*^) . (4.10) 

However, these results of averaging, change when we transform to spherical coordinates centred at O. When 
the components of the unit vector in Euclidean space are referred to the spherical coordinates, they are 
{v} ^v? ,n^) = (1,0,0) (with {x^ ^x^ ,x^) = (r, a, /3)), and then the averages are <ti'^vP>iT^= 5\5-[, <n'n^n'^>47r= 
5\S'!^5\. Explicit calculation done in Sec. 5 for the L T model will show that the general formulae (4.8) - (4.10), 
in particular (4.9), lead to incorrect results when applied at the centre of symmetry. 

Since (4.1) is supposed to follow from (4.7) by this kind of averaging, the argument above invalidates (4.1) 
at the centre of symmetry of the L-T model even as an approximate equation. 



5 The sign of qz' can be any 

We will now use (4.7) to estimate HS's deceleration parameter qy at the centre of the L T model without 
averaging over directions. Identifying the coefficient of Dj^ in (4.7) with the H^i of (3.3) we find from (3.3): 

^ The averages are defined by 
and similarly for <n^n^ n'°>4iT ■ 



47r Jo Jo 



on n-' sin a, 
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qy =1 + 



In the coordinates of (2.1) we have, with {x^,x'^,x^) = (r, cp): 



(5.1) 



(/ill, /i22, /i33) = ^Yli' ^^^^ ^ ' (^-^^ 

{Kii,K22,K33) = (^-f^^,-RR,t - RR,tsm^ i}^ , (5.3) 

^n|i = i + 2i; (^-^^ 

Kl2\2 = ^22|l = ^ = ^ = ^,t^,r - RR,tr, (5.5) 

sm w sm v 

the components not Usted are zero. 

At a general point of an L-T manifold, the unit spatial vector n' must obey hiju'^nP = 1, and so, in consequence 
of (5.2), its components are constrained by 



^1 __^ VTT2g 



At the centre of symmetry each vector n* is necessarily radial, so = = 0, and (5.6) gives 

n\0) = ^^l±^, (5.7) 

which is finite by (2.9) and (A.4). 

However, (5.6) and (5.7) imply that at the centre of symmetry the angle-average of (4.7) is discontinuous. 
Namely, as long as the vector n' is attached off the centre, its v} = vT component may point in the direction of 
increasing r (in which case we take (5.6) with the + sign) or toward decreasing r (in which case wo take — in 
(5.6)). Then the average comes out zero. But when the vector n* is attached at the centre of symmetry, 'n} may 
point only toward increasing r, so > on the whole sphere, while = n'^ = 0. Thus, (n^)"^ > everywhere on 
the sphere, while all other components oin^n-'n'' are zero. Consequently, ((n^)^)4^ > 0, while all the remaining 
averages are zero. Thus, in the spherical coordinates, {n^n-^n^)4^ is discontinuous at the centre of symmetry. 

Let us apply (5.1) at the centre of symmetry in the L-T model. Using (5.2) - (5.7) and 

Hs> = -Kijn'n^ lim (5.8) 

in (5.1) we find 

In order to clearly see the contribution of the inhomogeneous part of the geometry, we will separate R and E 
into the part that survives in the Priedmann limit and the part that disappears in that limit. Note that the whole 
calculation will bo exact, and our definitions for M, E and R will ensure they automatically obey the conditions 
of regularity at the centre (2.9). 

We define: 

R = r(S + V), limP = 0, (5.10) 
E = r'^{-k/2 + F), limF = 0, (5.11) 

1 — i-O 

where fc is a constant (the familiar FLRW curvature index, so it may be of any sign, also zero) and S{t), F{r) 
and 7-'(t,r) are functions {S being the FLRW scale factor). The terms rV and r-^F represent non-Friedmannian 
contributions to R and E, respectively; with V = F = the L-T model reduces to that of Friedmann. In 
calculating the limits at r — >^ we will assume that the derivatives of F and ts do not diverge too fast, that is 

Jim [r(^F,r,F,r^F,rr,tB^r,tB,rr)] = (0,0,0,0,0). (5.12) 

Since some of the functions will be complicated, we introduce the following abbreviations: 



y=Ji^-k + ^F, (5.13) 
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def 2{S+T)F,r _ 

k-2F 



k-2F 

(it can be verified that C = V,r), and the subscript will denote the limit at r 0, thus 



(5.14) 



Vo=J^-Mo_k, (5.15) 



(5.16) 



The formulae for the derivatives of R and their limits at r — > are given in Appendix A. We see from there 
that all the terms entering qy in (5.1) have nonzero and nondivergent limits at r — 0. 
For H^i we get from (5.8): 

lhnHy = ^. (5.17) 
Using (5.2) - (5.17) and (A.l) - (A.8) in (5.9) we now obtain: 

The terms proportional to F_r{0), f_B.r(0) and Fr(0)/fe can each have any sign, and so each one can make Q3'(0) 
negative. Note that F,r(0) (which comes from the non-Friedmannian contribution to the energy function) alone 
can cause 93' (0) < (when tB,r = 0); aiid so can ts^riP) alone (when F = 0). 

At the centre of spherical symmetry, q^i = = 54. This confirms that HS's result is not correct at r = 0, 
and a negative §4, as found by INN, is perfectly possible. 

The reason why HS concluded that 54 > is that they found {K^j^j^n^n-' n'')4^ = while averaging (5.1) which 
is not true at a spherical centre in spherical coordinates. In consequence, they dismissed all terms in (5.18) which 
come from K^j^/^n^n-' , i.e. all except the first one. With only the first term, one obviously gets limr-^o93' — 0- 

Since nothing can be concluded in general about the sign of qy (0), there would be no contradiction involved in 
certain L-T models imitating a negative FLRW deceleration parameter even if eq. (5.1) were exact. But we recall: 
this equation being approximate, using it to evaluate the sign of 93/ cannot lead to an unambiguous conclusion, 
so the alleged contradiction was a nonexistent problem. 

There is one more piece of evidence that HS's method of averaging is not universally correct. Suppose (5.1) 
is applied at the centre of symmetry of an L-T model. Then each component of the sum in (5.1) is invariant 
under the group 50(3). Thus, averaging over directions should be an identity operation and change nothing. 
This is indeed the case when one uses (5.3) - (5.5) for Kij and Kj^j^^., applied at r = 0, and (5.7) for n*(0); 
actually one gets (5.17) - (5.18) again. However, if one uses HS's prescription for averaging, then the whole term 
{K^j^l^n'^n■'n''){0) still drops out in spite of being spherically symmetric and nonzero initially. 



6 Problems with Flanagan's eq. (5) 

Flanagan's paper [18] addressed the question of whether superhorizon perturbations can cause apparent acceler- 
ation. He concluded that if their effect on q is negative, it must be too small to ho responsible for the acceleration 
of the universe. The paper was not intended to be applied to exact inhomogeneous cosmological models, so 
the author should not be blamed for the incorrect use made later of his equation. But the HS subcase of that 
approximate equation was used by VFW as if it were exact. 

In this section we follow F's reasoning in order to identify the approximations assumed along the way. 



6.1 Comments to F's reasoning up to eq. (13) 
The correct form of Flanagan's eq. (9) is: 

lf{x,x') = u°'{x) + u°''^{x)a.0{x,x') + ^u"^^{x)a.0{x,x')a--y{x,x')+O{s^). (6.1) 

The meaning of the symbols is as follows: x is the set of coordinates of the observer at the point P, x' is the 

set of coordinates of the light source at the point Q, u"(x) is the four- velocity of the cosmic medium at P 
(the observer is comoving), u°'{x,x') is the four- velocity of the medium parallel-transported from Q to P along 
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the light ray, a{x,x') is Synge's [21] world function for P and Q, and the coefficients are u°'^{x) = —u°'''^{x), 
M"/37(a;) = „";(/37)(a;). 

The quantity in (6.1) is apparently only needed to calculate the redshift, which F writes as 

1 + ^=^, (6.2) 

where ka is the tangent vector to the ray at P. Usually the numerator is taken simply at Q, without parallel- 
transporting it to P [7]. Nevertheless, it is easy to verify that (6.2) is equal to the usual formula. 

F then considers two families of light rays: one converging at P, with the afhne parameter s and tangent 
vectors k normalised so that kau" = —1 at P, and another one diverging from Q, with the afiine parameter A 
and tangent vectors I normalised so that lau" = — 1 at Q. In consequence of this, and of (6.2), the two afiine 
parameters are related by 

X = s{l + z). (6.3) 



6.2 The On approximation 

Below his (13), F writes: We choose the normalization of A so that A 1/X for A ^ near Q. This is an 
assumption equivalent to HS's (26) - our (4.2) {A is the same quantity as in (4.5)), and F makes use of it further 
on (see below our eq. (6.7)). This normalization also defines the units for A as the units of distance. 



6.3 F's equation (14) 

In the paragraph of F's paper that contains eqs. (14) - (16) the quantities referring to P and to Q are mixed 
up. The text below is our interpretation of this segment of F's reasoning, in which we take care about properly 
distinguishing between these two points. 

From the definition of the energy flux at P one obtains: 



'^^'^ =T^0{P)uf.u'^P = Ap^lau")l=Ap''{kc,u'')l/{l + zf = Ap^/il + zf, (6.4) 

\2 1 11 1 / 7 a \ 2 



dtd'^A 



whereby a rewriting error was corrected above: F's (kau") should read (lau") . The following information 
was fed into (6.4) at the consecutive equality signs: the definition of energy flux, the deflnition of the radiation 
energy-momentum tensor at Q, equation (6.2), the normalization of fe" at P. 

Integrating the energy flux at Q over a sphere of small radius Ag, one obtains for the luminosity at Q, again 
from the deflnition: 

^ = Aq^ (Za«")Q X 47rAQ' = ^ttAq^Xq''; (6.5) 

the simplification occurs because of the normalization assumed at Q. Putting (6.2) and (6.3) into the deflnition 
of the luminosity distance (which we denote here by Dj^ for consistency with the previous sections) one obtains 

^qAq = (6.6) 

In the next step F uses the assumed normalization of Aq to write the above as Dj^ = (1 -|- z)/A, but then 
gives up on it and uses (6.3) to rewrite (6.6) as follows: 

Dl^{1 + z^sqAq/Ap = (1 + zfsQ/{ApXQ) (6.7) 

(the second equality follows by the assumed normalization of A at Q). Then F cites Visser [22] for the result 

ApXq = 1 + 0(sq2). (6.8) 

Thus flnally, up to s^-terms: 

DL^{l + zfsQ, (6.9) 

which is the same equation as HS's (29). Just as in the HS paper, this equation is approximate, where in addition 
to HS's (26) (our (4.2)), Flanagan used one more approximation - our (6.8). 

Prom (6.9) and firom eqs. (1) and (11) in F's paper, F's equations (15) and (16) follow by simple substitution. 



6.4 Choice of averaging procedure 



Below his eq. (16), Flanagan notes that there is no unique way to choose the order in which the coefficients in 
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the Di, [z) equation should be averaged over directions, but dismisses this problem by saying: 

Thus the different averaging prescriptions give different answers. However the fractional differences are of order 
a /9 , which we have argued above is of order e and is small. 

This is correct in the context of his paper. However, this indicates that there is a further approximation involved 
in the reasoning, and this one is quite arbitrary and unpredictable. Namely, it depends on the will of the person 
doing the calculation exactly which quantity he or she wishes to average first. F himself uses Hq = (A~^) to 
obtain the -ffo of his eq. (4), but H^^ = (A) for computing his J. This is perhaps the strongest indication that it 
makes no sense to use the resulting final equation for determining the sign of any quantity. Note that the same 
problem exists for the HS derivation, but it was not mentioned in the HS paper. 



6.5 The problem with rotation 

The averaging over directions is defined only in a certain 3-space S3, not in spacetime. Without fixing S3, the 
directional angles of the light rays are not defined. When the rotation of the cosmic fluid is zero, as in HS's paper, 
S3 is the 3-space orthogonal to the fluid's flow lines. 

However, when the cosmic fluid has nonzero rotation, such S3 do not exist - not even locally, because the 
3-volume elements locally orthogonal to the flow lines, when followed around any flow-tube, refuse to connect up 
to a 3-space. One could consider a 3-space of constant time-coordinate, which is not orthogonal to the flow lines, 
but then the metric in this space is not the /i^^^ = stq,^ -|- MaW/3 assumed by F. Thus, in the presence of rotation, 
the whole calculation should be reformulated. Without that, for the averaging over angles considered by F there 
exists no space in which it occurs. 



7 VFW's weak singularity is not a singularity 

In the paragraph containing their eq. (2.18) VFW say: 
We expand the density (2.5) to second order in r as 

p{r,t) = po{t) + pi{t)r + P2 (i)r' + ©(r'^). (2.18) 

The weaJc singularity occurs when pi (t) is nonzero, in which case the gravitational held is singular since OR. — > ±00 
as r 0, where 71 is the Ricci scalar. In other words, second derivatives of the density diverge at the origin, 
independent of where observers may be located. This is true both in fiat spacetime and in the curved LTB metric 
when we have a density profile of the form (2.18). The singularity is weak according to the classification scheme 
of the literature on general relativity [28]. (Ref. [28] in this quotation is [6].) 

In truth, whether a curvature singularity is there is decided by curvature alone, and not by secondary con- 
structs like the d'Alembertian of the curvature scalar. We do not know about any physical or geometrical in- 
terpretation of this quantity, and VFW do not mention any, nor have they shown that their 'weak singularity' 
causes any problems at the origin. As we show below the curvature of the L-T model at the centre of symmetry 
is nonsingular provided that p is finite there. 

For the metric (2.1), the orthonormal tetrad components of the curvature tensor, in the basis defined by 

e° dt, R,rdr/^/T+2E, RM, =^ RsinMip, are 

2M M,r M 
Jtoioi = - B?Rr ' -"0202 = -Hoaos = 

i?1212 = i?1313 = --^ i?2323 = -— (7 1) 

«1212 «1313 ^3 «2323 ^3 • W-lj 

The quantities given above are all scalars, so any scalar polynomial in the curvature components will be a 
polynomial in the quantities given in (7.1). If these are nonsingular, then there will be no scalar polynomial 

curvature singularity. 

Let r = Tc be the radial coordinate corresponding to the centre of symmctrjr if = 0. The requirement of no 
point-mass at the centre implies M(rc) = 0. It can be seen that with M(rc) = R(t, re) = all we need to make 
the curvature nonsingular at r = re is a finite limit of M/R^ at r = re- We have, using (2.4) 

A^. f = }Tr. ^ = (4/3)-«^(^' -^)' (7-2) 

SO the curvature is finite aX r = Vc if and only if p{t, r) is. 
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As regards geodesic completeness, the geodesic equation is well behaved through the origin. Geodesies passing 
through the origin must be purely radial, and radial geodesies stay radial. The only Christoffel symbols that 
diverge at i? = are F^rO = r'^r<l> = R'/R, but these do not appear in the equations of radial geodesies, so 
cannot cause any problem. 

The condition OTZ ±00 is not a criterion for a singularity in any accepted sense. The paper of Tipler [6] 
contains a proposal of a definition of a strong singularity, but does not mention a 'weak singularity' anywhere and 
contains no 'classification scheme'. For a singularity to be 'weak' (i.e. not strong), it first has to be a singularity. 
An accurate description of the feature that VFW had in mind is 'the density profile is not through the origin'. 

Therefore, VFW's weak singularity is not a singularity, so its presence is no reason to dismiss models con- 
taining it. There is nothing wrong with a jump in the density gradient at the centre. Since our Galaxy is centrally 
concentrated and has a central black hole, a pointed density profile is a good approximation. The NFW density 
profile [23] for galaxy clusters is divergent at the centre. A jump in the density and its gradient exists, for example, 
at the surface of the Earth, thus we deal with a 'weak singularity' in everyday life. 



8 Transcritical solutions are generic 

A large part of VFW's paper is devoted to considering the inverse problem and arguing that only FLRW models 
have rays that cross the apparent horizon — what they call the critical point. We here dispel the impression they 
convey, that transcritical solutions are well-nigh impossible to find for non-homogeneous models that account for 
observations. We emphasise (a) that all L-T models with a decelerating expansion phase have apparent horizons 
that are crossed by large families of light rays, and (b) no cosmological model has been shown to solve the inverse 
problem (as distinct from the fitting problem) relative to real data. 

VFW correctly point out that the differential equations (DEs) of the inverse problem, i.e. the DEs determining 
the L-T model that reproduces given observational data, contain terms that would diverge unless a certain 
condition is met. They call such loci critical points, and they correctly claim that such points are generic to these 
DEs. This was actually pointed out by Mustapha, Hellaby and Ellis [25], solved by Lu and Hellaby [10], discussed 
as an observational feature in [11], and used to advantage by McClure and Hellaby [12]. These papers develop a 
numerical procedure for extracting the metric of the cosmos from observational data.^ 

In their analysis of the inverse problem, VFW consider only flat L-T models, which are a subset of measure 
zero in the family of all possible L-T models. They say that FLRW models provide examples of transcritical 
solutions (i.e. those that do not cause any divergence of the DEs), but they fail to find any other viable solution.^ 
They then conclude that physically reasonable solutions must be very exceptional indeed, and although they have 
not tried the E ^ Q case, they opine that generic solutions with E{r) 7^ would have all the singularities and 
observational problems they encountered, and in this case too transcritical solutions do not appear to be likely. 

In assessing the likelihood of finding a 'transcritical solution' it is important to understand the physics and 
geometry of the problem. For any given cosmological metric, the DEs that determine the path of the light ray and 
the variation of z and Dj^{z) along it (for example (2.5), (8.2), (2.6)) are free of critical points. However, the path 
followed by an incoming light ray depends on the geometry it passes through, and this affects the observations, 
z, Dj^{z) etc. All known expanding cosmologies with a deceleration phase have a past apparent horizon (AH~); 
relative to a given worldline (observer), it is the locus where incoming light rays overcome the cosmic expansion 
and start to make progress inwards. The maxirrmm in the area distance Dyim = ^m, occurs where an observer's 
PNC crosses her past apparent horizon. The fact that dR/dz = at this maximum is the cause of the critical 
behaviour. This behaviour is a generic feature of the L-T and Friedmann models. We express this in the following 
theorems. We consider only models that: (i) have a big bang and after 14 billion years are still expanding, so any 
recoUapse is to the future, and times well beyond the present are not considered; (ii) have no shell crossings [26]; 
(iii) are large, e.g. if the curvature is positive, then any spatial maximum where M,r = = R,r is well beyond 
observational range, which means M,r > and R,r > 0; (iv) if yl > its effect is only discernable in fairly distant 
observations. We will call such a model an rLTc: a realistic L-T cosmology. 

Theorem 1 In every rLTc, every light ray arriving at the origin is a transcritical ray. 



^ This problem persists even in the fluid-ray coordinates (observer coordinates); in [29], tlie differential equations (59)- 
(65) contain dC/dz in the denominator, and on the observer's past null cone w = wo, this is dR/dz here, so it is again 
zero at the apparent horizon. 

* An L-T model is fully specified by two physical functions. We agree that a random combination such as = and an 
arbitrarily chosen Dl{z), specified via rpRvviz) or V{z), may well produce an L-T model with unrealistic features. 
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Proof The areal radius R{t, r) obeys 



where I = +1 if i? is increasing (expansion) and —1 if it is decreasing (contraction). Worldlines with 2E > 
are ever-expanding and the rest arc rccoUapsing. Those that recoUapse reach their maximum R 
where 6M + 6ER + AR^ = and this maximum is < (3M/A)^/^. 

The incoming radial hght rays satisfy (2.5), so the variation of R along a ray is 



dR 
dr 



ar 



VTT2E I ' VlT2E dt ^ ' 



The past apparent horizon is the locus where dRjdr = for all rays, which implies 

6M = 3R- AR^ (8.3) 

and the past apparent horizon has £ = +1. For a given M > 0, (8.3) has either two, one or zero solutions in the 
range R > 0. The smaller R solution corresponds to the regular {A = 0) AH, and the larger one corresponds to 
the de Sitter horizon, but the two merge on the 3My/A = 1 worldline. 

Near the origin, R{t,0) = V t, regularity [30,28,27,16] requires R ~ M^/^ and 2E ~ M^/^, so by (8.1) 
R^t and by (8.2) dR/dr — > R^r there. This means 6M < 3R — AR^ for incoming rays near the origin. 
Following the rays outwards and back in time, (2.5) shows r and therefore M is increasing, and (8.2) shows R 
is increasing near the origin. Provided R.r stays positive as assumed, (2.5) may be integrated all the way to the 
big bang; and here we have M > and i? 0, so 6M > 3ii - AR^. Therefore the rays have crossed AH . This 
argument applies to every ray arriving at the origin before the big crunch, if there is a crunch. 

The DEs for the general L-T 'inverse problem' are given in [10,12], and it is evident that the critical points 
are where di?/dr = = dR/dz. If a ray crosses the past AH, then by definition it (a) has d7?/dr = there, and 
therefore (b) it is transcritical. □ 

If yl > then there could^ be worldlines with M > 1/{3\/A ). These worldlines never encounter a solution 
to (8.3), so every incoming ray crossing them has dR/dr < and d7i/dt > 0. In ever-expanding models, AH~ 
asymptotically approaches R = i/S/A (and M = 0) towards the future, and it separates incoming light rays that 
always have increasing R from those that reach the origin. Therefore there are rays emerging from the bang that 
neither cross AH~, nor reach the origin. 

Theorem 2 Every incoming radial light ray in every yl = L-T model with a big bang is a transcritical ray. 
Proof When yl = 0, (8.3) becomes the familiar 

R = 2M (8.4) 

and such apparent horizons have been well studied in [31,32]. In this case, every worldline encounters AH~, 
which is at -R > except at an origin, and every light ray starts on the big bang. Therefore, of the rays emerging 
from the bang and arriving at the crunch or i? = oo, only those starting at an origin may be said to have not 
crossed the past AH. Models usually have one origin, but closed models normally have two. The rays emerging 
from the local origin are outgoing, not incoming, and the rays emerging from an antipodal origin may eventually 
be locally incoming, but in an rLTc this would happen long after the present, if at all. □ 

Theorem 3 Transcritical rays of inhomogeneous models that reproduce the observed Da{z) do exist. 

Proof See e.g. [33,34]. The point here is that if a chosen model produces a Da{z) and D]^{z) that fits the 
observations for a given PNC, out past a maximum in D^, then that PNC is necessarily transcritical. □ 

Theorem 4 Any reasonable Dji{z) can be reproduced by a A = L-T model. 



^ In i? < models, it is possible tlie maximum M does not reacli 1/{3VA ). 
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Proof This is a special case of the theorem in [25] , which claims any pair of reasonable observational functions, 
Da{z) or Di^(z) and ij,n{z) can be reproduced by an L-T model. Here n{z) is the number density of sources 
(e.g. galaxies) in redshift space, and /x is the mean mass per source, which may also vary with z. Essentially the 
2 arbitrary functions in an L-T model allow one to reproduce 2 observational relations. In fact [11,10,12] have 
demonstrated how to solve the 'inverse problem' numerically for both observational functions, and shown that 
the conditions at the AH do impose a significant constraint that holds only at a single z value. But with only 
D[^{z) given, the 2 L-T functions allow the observations and the constraint to be fitted easily. □ 

Together these results show that there is a plentiful supply of transcritical solutions amongst the L-T models, 
and that there is a well defined procedure for extracting an L T model that fits the observed Da(z). 

More important, however, is the fact that VFW do not try to find an L-T model that reproduces real 
observational data, but the one that reproduces a chosen idealised function for Di^{z). 

VFW write as if FLRW models are known to solve the inverse problem. To date the 'inverse problem' has 
never even been attempted with real data. What is usually done is to solve the fitting problem: choose your 
favorite type of cosmological metric and find the free functions and parameters that make it best fit the data. 

Even models that arc arbitrarily close to reproducing the observations will fail this test. For example, take 
the 'observational' relations of a given FLRW model, add small random or systematic errors, and use this 'data' 
as input to the inverse problem. The DEs of the inverse problem will return the given FLRW model quite closely, 
but (8.3) will not be exactly satisfied where d7?/d2: = 0, so they will still diverge (see [12]). The chances of an 
FLRW light ray being exactly transcritical relative to real observational data are distinctly less than for an L-T 
model. 

The failure of a proposed solution to be transcritical does not indicate any physical problem with the under- 
lying L-T or FLRW model; at worst it indicates that a particular ray in a particular model does not (exactly) 
reproduce the given observational function(s). 

As discussed in [11,10,12,34,35], when attempting to reproduce two observational functions with an L-T 
model, the AH constraints do provide a significant challenge. In overcoming this challenge it was shown [11,12] 
the AH can give us useful observational information. 

9 Conclusions 

We have corrected a number of inaccurate statements appearing in the HS [2], Flanagan [18] and VFW [1] papers. 
At least two of these misunderstandings are widespread and therefore it is worth putting the record straight. 

The first error-correction to be made is that a conical density profile at the origin is not a singularity in any 
accepted sense, and there arc no physical problems associated with it. In some situations a conical profile may 
be a good model of structures that are strongly centrally concentrated. 

Secondly, critical points in the inverse problem are generic. When solving the 'inverse problem' as defined in 
Sec. 8 for a realistic L-T or FLRW model, and finding that the differential equations diverge at some redshift, 
one must be aware that this is an inevitable consequence of the apparent horizon, where the diameter distance 
is maximum. The chances that real observational data exactly satisfy the full transcriticality conditions for any 
given model, homogeneous or not, are essentially zero. Such divergent behaviour does not indicate any physical 
problem with the underlying L-T or FLRW model, and transcritical solutions in a variety of models that closely 
reproduce the observed D]^{z) do exist. Actually this point seems not to be well understood and other examples 
of this misunderstanding can be found in the literature (see, e.g. [3,36]). 

Our point three is that one should be careful when extending FLRW definitions of some functions to other 
cosmological models. In most such cases, when inhomogeneities are present, the FLRW parameters do not mean 
the same thing as in the homogeneous case. 

It is very common to get acceleration and supernova dimming mixed up. The interpretation of the supernova 
dimming as an accelerated expansion (which is due to the use of FLRW models) has firmly taken root, which is 
why some have focussed on acceleration in the fluid expansion. In the L-T models, the 'acceleration' defined from 
z{D[^) using (3.3) exhibits a behaviour different from the qi of the fluid flow defined by (2.1), even though they 
coincide for FLRW models. We hope our paper will help establish the difference between 'accelerated expansion' 
and 'supernova dimming' in the observed magnitude-redshift relation. 

We summarise and discuss different misunderstandings related to inhomogeneous cosmology in Table 1. 
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Table 1 Suiiniiar\" of iiiist'ouccpl iuiis aud iiiisuij(l(*rs(an(liii.i;s r('lat(Ml U) iulioiiioyciicous t*(jsiii(jl()t;\". 



Misconceptions and misunderstandings 



Corrections and clarifications 



Weak singularity 



DTZ — >■ oo is a singularity. 

Models with \U]Z = oo are unphysical. 



OTZ — >• oo is not a curvature singularity and has no physical 
interpretation. Many objects have = oo, among them the Earth at 
its surfa<;e. 



There are general theorems that prohibit go < where 
go = q{z = 0) presented in Refs. [18,2]. 



There is a local singularity in models with go < 0. 



Deceleration parameter 

There are two distinctive deceleration parameters: g"*"* based on a 
Taylor expansion of the luminosity distance and the invariantly defined 
q'"" which measures the acceleration of expansion. If A = = u)„i, = 
and p + 3p > then g'"" > 0, however for the same case g"*"* may be 
negative. F's & HS's equations relating g"*"* to the flow invariants are 
approximate and coordinate-dependent, thus they do not exist as 
covariant laws. The approximately correct intermediate relation does not 
exclude q°^^ < in the L-T model. 

DTZ — > oo when go < and the observer is at the centre of spherical 
symmetry. However, this is not a singularity. Away from the the centre 
of spherical symmetry go < does not imply divergence of DT?.. 



Other singularities arise in models with q{z) < 0. In such 
models there is a class in which 



l+z dH{z) 



H{z) dz 



q{z) = -l 



oo at some z, where 



g defined in this way diverges when ^ (^y+z) ~ ^ ^^^^ ^® '^'^^ ^ 
singularity. In fact, if one applies these definitions to a zero- /I dust 
FLRW model, one obtains regions where c; < when Qf^Q > 0.6. 
This is because such H{z) only applies when A: = 0. 



In solving for the model that gives a selected Dij{z), one 
encounters a 'pathology' or 'critical point', beyond z ~ 1 
and the solution generally breaks down there. 



Only FLRW models have "transcritical" solutions. 



Inverse problem 

The 'critical point' is the apparent horizon, where the diameter distance 
is maximum. It is a general property of expanding decelerating /I = 
models such as L-T and FLRW, long known in the FLRW case. This 
point requires special treatment, but has useful properties [11,10,12]. 

"Transcritical" rays, that cross the apparent horizon, are generic in L-T 
and FLRW models. With small errors in observational data, all L-T and 
FLRW models will technically fail the apparent horizon conditions, but 
this can be fixed [12]. Arbitrary choices such as E = and a given 
Di^ {z) may well fail to be transcritical. 



Inhomogeneous models are exposed to singularities. 
Non-singular models are exceptional and rigid, hence 
unable to account for observations. 



gg*"^ < is essential to account for observations. 



Fitting observations 

Inhomogeneous models like the L-T models include the FLRW models 
as a subcase. Thus, if the FLRW models are considered good enough for 
cosmology, then the L-T models can only be better: they constitute an 
exact perturbation of the FLRW background, and can reproduce the 
latter as a limit with an arbitrary precision. 

Derivation of g"*"" in an inhomogeneous model involves approximations - 
one of them linearity, so the sign of g"*"* cannot be determined. One can 
have a very good fit to observations with q°^^ > [38,39]. 



A Behaviour of derivatives of R{t, r) in the neighbourhood of the symmetry centre 

The formulae below are obtained using the coordinate r introduced in (2.8), the symbols introduced in (5.10) - (5.11), and 
eq. (2.2). We show the intermediate expressions in order to demonstrate that all terms that could cause divergencies cancel 
out before the limit is taJjen. The first two equations follow trivially from (2.2). 

l2M 

= V:^ + 2£; = ry-^^0, (A.l) 
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R.tt = -- 



M 



Mor 



0, 



R2 (5 + ^)2 r^-O 

To evaluate the behaviour of -R^^ and R^rr in the neighbourhood of r = we use the following equation [37]: 



R., 



\ ^ r/3 E^r M,r \ , 



R.i 



(A.2) 



(A.3) 



The above applies also with S = if we neglect the E^r/E terms (as verified by direct calculation). Prom (A.3) we obtain: 



S + V + rC 



(A.4) 



The next equations below follow by consecutively differentiating (A.l) - (A.3) by r and using (A.3) - (A.4) in the result: 



R,ttr = — 



M, 



2MR, 



Mo 



R.i 



R,t \ R 



Mr MR, 



+ 



R^ 



+ E, 







2MorC Mo 
MqC 



V - 



V 



Vo, 



(A.5) 
(A.6) 



1 /Mr MRr ^ V 1 / 



Mrr 2MrRr MR, 



Mo 



V I {S + V) 
Mo 



1 



r-!-0 Vo I 52 



-R,rr+2Fr+r 
-R,rr-(0) +2F,T.(0) 



R R2 

2Mo/:2 1 



il2 



■ + ■ 



2MR,, 

'rF 



+ E, 



(s + vy 



r 

V3 



MoC 



(5 + ^)2 



(A.7) 



R 1 



1 



-k/2 + F 



-iS + V) + -V{t-tB) 



2ts,r + r tB.rr + 



F.rtf 



2 -k/2 + F 

2F,r(0) 



2(-fc/2 + F) 



{t-tB)-tB,, 



2Fr + r Frr - 2 



MoC 



-k/2 + F 



•' (5 + ^)2 
[-25 + 3 (t - ts(0)) Vo] - 2ts,^(0)yo- 



(A.8) 
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